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Recently published measurements of the branching ratios B{ip{lS) — > 7770(15')) and B{ip{2S) — > 
7 r?c(lS)) by the CLEO collaboration are examined in the context of a potential model that includes 
both relativistic and one-loop QCD corrections to the quark-antiquark interaction. The prediction 
for the width V{tl){lS) — > 7770(15*)) is in excellent agreement with the new data but the prediction 
for r(7/)(25) — > 7770(15)) is too small. In an effort to understand this discrepancy, we derive an 
upper bound on r(?/>(25) — > 7770(15)) and point out its experimental value saturates this bound. 

PACS numbers: 12.39.Pn 



1. INTRODUCTION 

In a recent publication [l| , the CLEO collaboration reports new measurements of the branching ratios of 
the charmonium radiative decays ^{IS) — > jric{lS) and 'ip{2S) jric{lS). The new values are B{ip{lS) — > 
7770(15)) = (1.98±0.09±0.30) X 10-2 and S(^/'(25) ^7770(15')) = (4.32±0.16±0.60) x IQ-^. These results 
imply radiative widths ofr^^p{i:{lS) -^j?]c{lS)) = 1.85 ±0.28 and rexp(V'(25') -> 7 ?7c(15')) = 1.41 ±0.21, 
both of which are larger than the current Particle Data Group values [3] • 

Magnetic dipole (Ml) transitions of this type particularly sensitive to the details of the charmonium 
radial wave functions and, as such, these data provide an important check on model calculations of radiative 
transitions. There are numerous approaches to these calculations including lattice QCD heavy quark 
effective theory Q , inclusion of hadronic loop effects @ and potential models 0] . In this note^we describe 
the results of comparing the new data with the potential model calculations contained in Ref. [y] . 



2. Ml TRANSITIONS 

In the dipole approximation, the width for the radiative transition Tp(nS) — > jr]c{n'S) is given by 



3 mj. M^,(„s) 



where {n'is'\n£s) denotes the radial integral 

{n'is'\nis)= I drr'Rr,is'{r)RnUr), (2) 
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where lu is the photon energy and -E^cCn'S) is the energy of the recoiling rjc- Here s denotes the initial spin 
(s = 0, 1), s' the final spin and s' = s ± 1. In a model in which the radial wave functions used to compute 
the Ml matrix elements are obtained using a Hamiltonian that does not contain spin-dependent terms, the 
£ = singlet and triplet radial wave functions corresponding to different radial excitations are themselves 
orthogonal. However, in general, £ = singlet states (s = 0) and £ = triplet states (s = 1) are orthogonal 
by virtue of their spin wave functions so there is no reason why the singlet radial functions i?n'oo('') should 
be orthogonal to the triplet radial wave functions i?„oi(f) when n' ^ n. Given this, the radial wave function 
of the r]c{lS) obtained in our non-perturbative treatment can be expanded in terms of the ip{nS) radial 
wave functions wave functions as 

OG 

n=l 

with the C„'s given by 

/■oo 

C„ = / drr2i?„oiW-Rioo(r-). (4) 
Hence, the amplitudes (100|n01) are just the C„'s and, since the rjci^S) radial wave function is normalized, 

oc 

Y.\Cn\' = l. (5) 

Given a model that adequately describes r{tlj{lS) — > jrjcilS)), the value of [Cip can be used to obtain a 
bound on r{ip{2S) — > jr]c{lS)) by noting that Eq. ([5]) can be written 

oo 

i-i^ii'==Ei^"i'^i^2r- (6) 

n=2 



To see if the bound in Eq. ^ is particularly restrictive, we use the results of the potential model described 
in detail in Ref. Briefly, the model uses a Hamiltonian consisting of a short distance potential that 
includes all j <? relativistic corrections and one-loop QCD corrections, a linear confining potential that is 
a mixture of scalar and vector contributions together with their j <? corrections and a relativistic kinetic 
energy term. The charmonium spectrum is obtained using a variational technique that provides explicit 
forms for the radial wave functions. The amplitudes in Eq. Q can then be evaluated and the predictions for 
the radiative widths calculated using Eq. ([1]). 

The radiative widths calculated this way using the non-perturbative approach of Ref. @ are 

r(V'(15) ^7^c(15)) = 1.84keV, (7) 
r(Vj(2S') 7?7c(15)) = 0.44keV. (8) 

From Eq. ([7]), the width of the V'(l'S') — > 7 77c(15') transition is well described by the model, but Eq. ([5]) 
shows a discrepancy of about a factor of 3. Faced with this discrepancy, one might ask if it is possible to 
modify the 'ip{2S) — > 777c (IS') amplitude enough to achieve agreement with experiment without violating 
probability conservation. The extent to which this can be done is controlled by Eq. The calculated value 
of |Ci|2 leading to the resuh in Eq. © is |Ci|^ = 0.9958, so Eq. © becomes 

< 4.196 X 10-^ (9) 

whereas the calculation leading to Eq. ([5]) gives |C2|^ ~ 1.571 x 10""^. Rescaling the result in Eq. ([5]) by the 
ratio then implies 



T{'iP{2S) ^jricilS)) < 1.20 kcV 



(10) 
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Interestingly, the experimental value. Tcxp{ip{2S) — » 7?7c(lS')) = 1.41 ± 0.21, saturates this bound within 
errors. Before drawing any conclusions, it must be remembered that the bound is obtained using the dipole 
approximation, which amounts to replacing the complete Ml amplitude, 

/•oo 

A{Ml)=j drr^R„,es'{r)Rnes{r)MY)^ (H) 

with Eq. ([2]) by using jo{x) — > 1. This is justified as long as the correction from the next term in the 
expansion of jo (3-), a;^r^/4!, is small. In the present case, the wave functions used in the variational 
calculation fall off exponentially at large r with a scale parameter i? = 1 Gev~^. Using R to estimate 
the size of the correction, we find 5.32 x 10^^ for the ip{lS) — s- 777c(l'S') transition and 1.69 x 10~^ for 
ip{2S) — > 7 77(7(15) transition. This indicates that deviations from the dipole approximation do not materially 
affect the validity of the bound on r{ip{2S) 777,7(15') given in Eq. PH)) . 



3. CONCLUSION 

Within the confines of the dipole approximation, a bound of the type in Eq. ^ on the Ml transition 
amplitude from the first excited triplet state to the singlet ground state can always be obtained. Furthermore, 
the transition amplitudes in any sensible model calculation will satisfy this bound. 

So, is there reason to be concerned that Tcxp{tl){2S) — > 777,7(15')) just barely satisfies the bound associated 
with a particular model that accurately predicts Tcxpiipi^S) — > "f ric{lS))7 From the potential model point 
of view, the value of rexp('0(25) — > 777^(15)) is uncomfortably large because all the ingredients needed 
to evaluate Eq. ([T]) are tightly constrained by the requirement that the variational calculation accurately 
reproduce the observed charmonium spectrum Q. Once this is accomplished, the |C„p's are determined 
and those with ji > 2 are not zero. The challenge posed by the new data is one of seeing whether it is possible 
to refine the wave functions in such a way that IC2P can be increased within the constraints of probability 
conservation without sacrificing the quality of the overall fit to the charmonium spectrum. 
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